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Abstract 

^ ■ We develop a perturbative QCD factorization theorem which is compat- 

5h ' ible with effective field theory. The factorization involves three scales: an 

infrared cutoff of order Aqcd, a hard scale of order the B meson mass, and 
an ultraviolet cutoff of order the W boson mass. Our approach is renor- 
malization group invariant and moderates the scale-dependence problem in 
effective field theory. Applying this formalism to exclusive nonleptonic B 
meson decays, we clarify the controversy over the Bauer-Stech-Wirbel pa- 
rameters 02/0.1 for charm and bottom decays. It is found that the nonfac- 
torizable contribution plays an important role in the explanation of the sign 
and magnitude of 02/01. 
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Nonleptonic heavy meson decays are difficult to analyze due to the compli- 
cated QCD corrections. While the semileptonic decays involve only conserved 
currents, the nonleptonic decays are described by four-quark current-current 
operators, which is part of the low-energy effective Hamiltonians for the W 
boson exchange. For example, the relevant operator for the B Dtt decays 
is 

AGf 

H = -^V,tV:,{cLl,h){dLYuL) , (1) 

The QCD corrections will generate operator mixing among these operators. 
The mixing, characterized by Wilson coefficients, depends on an arbitrary 
renormalization scale /x. The full effective Hamiltonian related to Eq. ([^) 
can be written as 

HeE = -^V,bV:Aci{fi)0, + CMO2] , (2) 

with 

Oi = (cLlMidLr^L) , O2 = {dLl^,hL){cLYuL) . (3) 

In Eq. ci and C2 are the Wilson coefficients, whose evolution from the 
W boson mass Mw down to a lower scale is determined by renormalization- 
group running [||. Though the Wilson coefficients are /i dependent, physical 
quantities such as decay amplitudes are not. In principle, the matrix elements 
of the four-fermion operators contain a /i dependence, which exactly cancels 
that of the Wilson coefficients. In practical applications, however, various 
schemes are needed to estimate the hadronic matrix elements, and the esti- 
mates are usually yU independnt. Hence, the decay amplitudes turn out to be 
scale dependent. Take exclusive nonleptonic B meson decays as an example, 
to which the conventional approach is the Bauer-Stech-Wirbel (BSW) fac- 
torization approximation 0. It is assumed that nonleptonic matrix elements 
can be factorized into two matrix elements of (axial) vector currents. Since 
the currents are conserved, the matrix elements have no anomalous scale de- 
pendence. Presumably /i should be set to the dominant scale of the matrix 
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elements. However, the matrix elements involve both the heavy quark scale 
and the small hadronic scale. Naively setting fi to the heavy quark mass 
will lose large logarithms associated with the hadronic scale. It is then quite 
natural that theoretical predictions are sensitive to the scale we choose P, |]. 

To circumvent this problem, a phenomenological approach is adopted to 
bypass the strong scale dependence. The Wilson coefficients Cj are regarded 
as free parameters, and are determined by experimental data In this 
model, two equivalent parameters ai = Ci+C2/Nc and a2 = C2+ci/Nc describe 
the external and internal VT-emission amplitudes, respectively. However, the 
evaluation of the hadronic form factors usually involve some ansatz and 
thus the extraction of ai and 02 is model dependent. It is also found that 
only when the experimental errors are expanded greatly does an allowed 
domain (oi, 02) exist for the three classes of decays B'^ — ^ D^*^^, /)(*)o 
and B- -> [|. In addition, a negative 02/01 and a positive 02/0.1 are 

concluded from the data of charm and bottom decays [Q, [7|, respectively. 

It was shown recently that the perturbative QCD (PQCD) approach is 
applicable to heavy meson decays p, |^, to which the heavy mass provides 
the large momentum transfer. The breakthrough is due to the all-order 
Sudakov resummation of large radiative corrections, which suppress contri- 
butions from the long-distance region. This formalism, taking into account 
the evolution from the typical scale of the hard subprocesses characterized by 
the heavy meson mass to a lower hadronic scale, is /x independent for semilep- 
tonic decays. In this letter we shall develop a PQCD formalism based on the 
effective Hamiltonian in Eq. (^, which further incorporates the evolution 
from Mw down to the hard scale. This three-scale factorization theorem, 
being fj, independent, does not suffer the scale-setting ambiguity mentioned 
above. We apply this formalism to two-body nonleptonic B meson decays 
such as B ^ Dtt. Without any free parameter, our prediction agrees well 
with experimental data. 

We ffist illustrate the main idea of PQCD factorization theorems by con- 
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sidering one-loop QCD corrections to a generic decay process through a cur- 
rent. These corrections are ultraviolet finite, since the conserved current is 
not renormalized. However, they also give rise to infrared divergences, when 
the gluons are soft or collinear to light partons. The factorization is imple- 
mented to isolate these infrared divergences associated with the long-distance 
physics. 

Radiative corrections that produce infrared divergences are classified into 
the reducible and irreducible types []T0|. Irreducible corrections contain only 
single soft logarithms and is absorbed into a soft function U, which cor- 
responds to the nonfactorizable soft corrections in B meson decays in the 
literature 0. These corrections cancel asymptotically |jlO[, and thus are ex- 



pected to be small. They will be neglected here {ie U = 1) and studied in a 
forthcoming work. Reducible corrections contain double logarithms from the 
combination of soft and collinear divergences, which can be absorbed into a 
wave function (j){P,b,fi) and explicitly resummed into a Sudakov factor [|T^], 

0(P,6,/i)=exp[-s(P,6)]0(6,/i). (4) 

b is the conjugate variable of the transverse momentum, which will be ex- 
plained later, and 1/6 can be regarded as an infrared cutoff. 

To factorize a one-loop correction, we divide it into two terms as shown 
in Fig. 1(a). The first term, with eikonal approximation for fermion prop- 
agators, picks up the infrared structure of the full diagram. Being infrared 
sensitive, it is absorbed into U or (p, depending on which type the one-loop 
correction is. The second term, with a soft subtraction, is infrared safe. It 
has the same ultraviolet structure as the full diagram and can be absorbed 
into a hard scattering amplitude H{t, /i), where t denotes the typical scale of 
the hard decay process. We then get the 0{as) factorization formula shown 
in Fig. 1(b), where the diagrams in the first parentheses contribute to H. 

The presence of fi implies that both and H need renormalization. Let 'j^ 
be the anomalous dimension of d). Then the anomalous dimension of H must 
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be — 7(^, because the full diagram does not contain ultraviolet divergences. 
Their fi dependence can be calculated by RG, 



0(6, /i) 



(j){b, 1/6) exp 
H{t,t) exp 



— 70 «s /i 

1/6 /i 



7</>(a«(/^)) 



(5) 
(6) 



Equation (|[) describes the evolution of from 1/6 to an arbitrary scale fi, 
and (^ describes the evolution of H from fi to t. The contribution character- 
ized by momenta smaller than 1/6, ^e., the infrared divergence, is absorbed 
into the initial condition 0(6, 1/6), which is of nonperturbative origin. The 
convolution of H with is then /x independent as indicated by 



H{t, fj.)^{b, /i) = H{t, t)(j){b, 1/6) exp 



— 7<^(«s(/i)) 
i/b /i 



(7) 



In this way all the large single logarithms are collected in the exponential. 

Indeed the effective Hamiltonian in Eq. (Q) can be constructed in a sim- 
ilar way. Consider now a typical one-loop QCD correction to the W boson 
exchange diagram Fig. 1(c). We express the full diagram, which is ultravi- 
olet finite, into two terms as shown in Fig. 1(c). The first term, obtained 
by shrinking the W boson line into a vertex, corresponds to the local four- 
fermion operators Oj. It is absorbed into a hard scattering amplitude H(t, /i), 
with a typical scale t <^ Mw, since gluons involved in this term do not "see" 
the W boson. Its dependence on Mw is limited to the 1/M^ factor of the 
four-fermion operators. The second term, characterized by momenta of or- 
der M\Y, is absorbed into a "harder" function Hr{Mw, jj) (not a scattering 
amplitude), in which gluons do "see" the W boson. 

We obtain the 0{as) factorization formula shown in Fig. 1(d), where 
the diagrams in the first parentheses contribute to iJ^j and those in the 
second parentheses to H. Note that this formula in fact represents a matrix 
relation because of the mixing between operators Oi and O2. Solving their 
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RG equations, we derive 
Hr{Mw, ^i)H{t, jj) = Hr{Mw, Mw)H{t, t) exp 



— iHriasifi)) 



(8) 



where the anomalous dimension of is also a matrix. We emphasize 
that the factorization in Eq. is not complete because of the presence of 
infrared divergences in H. The exponential can be easily identified as the 
Wilson coefficient, implying that fi in c{fi) should be set to the hard scale t. 
Without large logarithms, Hr{Mw, Mw) can now be safely approximated by 
its lowest-order expression H^'^^ = 1. 

We are now ready to contruct a three-scale factorization theorem by com- 
bining Eqs. (0) and (P). Consider the decay amplitude up to 0{as) without 
integrating out the W boson. We first factorize out the infrared sensitive 
wave functions as described above. Though devoid of infrared divergences 
(the nonfactorizable soft corrections have been neglectd here), the hard part 
still invloves two scales t and M^r. The factorization in Fig. 1(d) is then 
employed to separate these two scales, and Hr can be moved out of the hard 
part, a step valid up to 0{as)- We identify the remaining diagrams, includ- 
ing the four-fermion amplitude and the associated soft subtraction, as the 
hard scattering amplitude H, since it is free of infrared divergences. The 
anomalous dimension of H is given by •yn = —{.l<t> + iHr)- We thus get the 
three-scale factorization formula 



Hr{Mw, f^)H{t, fj.)(j){b, fx) = c{t)H{t, t)(f){b, 1/6) exp 



, — 7<^(a.(/^)) 

1/6 /i 



(9) 



with the Wilson coefficient c(t) given by the exponential factor in Eq. (P). 
The two-stage evolutions from 1/6 to t and from t to M\y are both included, 
and the final expression is fi independent. 

The above conclusion is quite natural from the effective field theory ap- 
proach ||ll|. An effective field theory is constructed for a scale fi < M^y by 
integrating out the W boson at /i = Mw Matching corrections are deter- 
mined by the matching condition requiring that the low-energy light-particle 
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Green functions of the two theories be equal. The effective theory is then 
evolved by RG running from fi = M^y to a lower scale, which insures that the 
amplitudes are fi independent. The scale /i in a continuum effective field the- 
ory is actually a scale to separate the long-distance from the short-distance 
physics with the physics above the scale fi absorbed into the coefficients in 
the effective Hamiltonian, such as the Wilson coefficients Ci_2(/i). This idea 
is identical to the PQCD factorization theorem. The effective field theory 
constructed this way has exactly the same low-energy behaviour as the full 
theory, including infrared divergences, physical cuts, and etc. Thus the in- 
frared divergences in the decay amplitudes calculated using the effective field 
theory can be factorized in the same way as the full theory. The factorization 
formula for the /i independent amplitude is identical to Eq. (H), 

c{Mw,f^)H{t,fi)<P{b,fi) , (10) 

with the Wilson coefficient c identified as Hr- 

We now apply the above formalism to the nonleptonic decays B{Pi) — >■ 
D{P2) T^{P^)- The decay rate can be written as 

r = -^Gl\VMVu,?Ml^-^^^\M\' , (11) 

with r = Md/Mbi Mb (Md) being the B (D) meson mass. In the rest 
frame of the B meson. Pi has the components Pi = (Mb/v^)(1, 1, Ot). 
The nonvanishing components of P2 and P3 are respectively P<^ = Mb/V^, 
P2 = rMD/V2, P3+ = 0, and P3 = (1 - r^)MB/V2. Let ki{k2) be the 
momentum of the light valence quark in the B (D) meson and ^3 be the 
momentum of a valence quark in the pion. These /c's may be off-shell by 
the amount of their transverse components kx of order Aqcd- We define the 
momentum fractions x as Xi = fcf /Pi~, X2 = k2 / P2 ■, and 0:3 = k^ / P^ . The 
transverse momenta kiT play the role of an infrared cutoff in our analysis. 
To leading power in 1/Mb, the factorization formula for Ai in the trans- 
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verse configuration space |T0[ is written as 



[dx] / [d^b]0ij(a;i,6i,l/6i)0D(a;2,62,l/&2)07r(a;3,&3,l/&3) 



M = 

Jo ' ' Jo ' ' ' . , . . , . , . 

xc{t)H{xi, k, t) exp[-5(a;„ k)] (12) 

with [dx] = dxidx2dx3 and [d'^h] = d'^hid'^h2d'^h3. The Sudakov factor e^"^ is 
the product of in Eq. (^) and the exponential in Eq. @ from each wave 
function. In the analysis below we shall neglect the b dependence of the wave 
functions 

Without large logarithms, the hard part H can be reliably treated by 
perturbation theory. To leading order in a^, the hard part for the decay 
B~ — ^ D^Ti~ consists of four sets of diagrams shown in Fig. 2. The diagrams 
in Fig. 2(a) correspond to the external W emission p, while those in 
Fig. 2(b) to the internal W emission. They have been calculated using the 
PQCD formalism in ^ without including the Wilson coefficients. Denote 
their contributions to the amplitude as and A^5. It is easy to find 
that the Wilson coefficients associated with M.^ and M.h are respectively ai 
and a2- Readers are refered to [^J for the complete formulas of Ada and M-f,. 

Diagrams in Fig. 2(c) and 2(d) are absent in the factorization approxi- 
mation and will be called the nonfactorizable diagrams. Fig. 2(c) leads to 
the amplitude M.c 

Mc = 32Jm,7rCFV^MlGF f\dx] r hdhb2db2^B{xi)(l)Dix2)Mx3) 
* JO JO 



X 



«s(ti)^e-^-(-'^^)(xi-X2-X3(l 



r'))h^^\x,,k 



+«s(t2)^e-^-(^-''')(l - (xi + X2){1 - r'))hf\x,, h) 



(13) 



The functions W\ 



1 and 2, are given by 



h^p = [e(b^-b2)K^{AMBbi)h{AMBb2) + d{b2-bi)Ko{AMBb2)Io{AMBbi)] 
' K^{B,MBb2) 



X 



T^o 



\BAMBb2 



for Bj > 
for Bj < 



(14) 
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with = XiXsi^l — r^), Bf = (xi + X2)r^ — (1 — Xi — X2)xs{l — r^), and 
B2 = {xi — a;2)x3(l — r^). The Sudakov exponent Scj is written as 

Scj = s{XiP^, hi) + S{X2P^, h) + - X2)P2^, 62) + S{X^P^, 63) 

+.((l-.3)P3-.M-^i:ini^, (15) 

with 63 = 62) A = (33 — 2r2j)/12 and rij = 4 the number of flavors. The scale 
tj is chosen as tj = max^AMs, \Bj\MB, I/&1, 1/&2), and A = Aqcd is set to 
0.2 GeV. The Sudakov suppression described by the factor exp(— ^c) then 
warrants that main contributions come from the small b, or large t, region, 
in which as{t) is small, and thus the perturbative treatment of the hard 
part is reliable. The amplitude A^^ is obtained from Fig. 2(d) accordingly. 
The amplitudes for the decay B^ D^tt~ can be derived in a similar way. 
However, it is found that only the external VT-emission contribution, the 
same as Ada, is important. Therefore, we shall not give its expression here. 
The wave functions are chosen as 

M-) = ^Ml-x){l-2x) , M^) = J^Ml^ + Csn{l-x)^ 

(16) 

where = 132 MeV is the pion decay constant. Nb = 650.212 and Cb = 
—27.1051 correspond to the B meson decay constant /b = 200 MeV. No is 
determined by the D meson decay constant fo = 220 MeV, and is fixed 
by data for the decay B^ D^n^ All other parameters are referred to 

i- 

The experimental data of the branching ratios are Bq = B{B^ D^ti~) = 
(3.08 ± 0.85) X 10-3 and i3_ = B{B- D^n-) = (5.34 ± 1.05) x 10"^ [|l|]. 
Our predictions using the original Hamiltonian in Eq. (|l|), ie Ci = 1 and 
C2 = 0, are Bq = 3.08 x 10"^ and B- = 5.10 x 10"^. If the three-scale factor- 
ization formulas based on the effective Hamiltonian is employed, we obtain 
Bq = 3.08 X 10"^ and B- = 5.00 x 10^^, which differ from the previous results 
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only by 2%. The slight variation justifies the analysis of exclusive nonlep- 
tonic B meson decays performed in using Eq. A careful observation 
reveals that when the evolution of the Wilson coefficients is taken into ac- 
count, the amplitude J^b, which is proportional to 02 (t), becomes smaller, 
while Aic, which is proportional to Ci{t)/Nc, becomes larger. Hence, the 
two changes cancel each other, and the total decay rate remains almost the 
same. Aid is less important because of the pair cancellation between the two 
diagrams in Fig. 2(d). Our calculation also indicates that the nonfactoriz- 
able contribution Ai^ is substantial, and the limit of the BSW factorization 
approximation. That is why the naive choice of 01^2 = 01,2(^5) in the BSW 
model fails to explain the data. 

Applying the three-scale factorization theorem to the mode K^n~ , 

we obtain the predictions i3(Z)° K+tt') = 4.05% and B{D- -> K^tt') = 
2.67%, consistent with the data (4.01 ± 0.14)% and (2.74 ± 0.29)%, respec- 
tively. With the running scale t reaching below the c quark mass. Ait be- 
comes more negative and overcomes the positive contribution of Aic- This 
explains the observed destructive interference of the external and internal 
W-meson emission contributions absent in the B meson decays. Hence, the 
nonfactorizable diagrams play an important role in the explanation of the 
charm decay data. It is clear that a PQCD formalism based on the origi- 
nal Hamiltonian without Wilson coefficients [^, |l^ can not account for this 
change of sign in the charm decays. From the above analysis, we suggest 
that ci,2 could be regarded as Wilson coefficients as they originally are, if the 
scale is chosen properly, instead of as fitting parameters in the BSW model. 
That is, the controversy over the extraction of 02/(21 from the bottom and 
charm decays does not exist in our theory. 

The scale dependence of our formalism can be tested by substituting 2t 
for t in the factorization formula. It is found that the prediction decreases 
by only 5%. In the conventional effective field theory, the substitution of 
by 2Mb for the argument of the Wilson coefficients ci,2 results in a 10% to 
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20% difference 0]. Hence, the scale-setting ambiguity is moderated in the 
three-scale factorization theorems. In conclusion, our formalism provides a 
more sophisticated choice of the scale, which takes into account the additional 
low-energy dynamics in the mesons. This approach is expected to give more 
definitive predictions, when it is applied to inclusive nonleptonic B meson 
decays. This subject will be discussed elsewhere. 

Our formalism can also be applied to the decays B J/'^K^*\ With 
the three-scale factorization formulas and the inclusion of nonfactorizable 
contributions, both the decay rates and the fraction of the longitudinal mode 
can be explained. The details will be published in a separate work. 

This work was supported by the National Science Council of ROC under 
Grant No. NSC-85-2112-M-194-009 (for H. L.) and NSC-85-2112-M007-029, 
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Figure Captions 



Fig. 1. (a) Separation of infrared and hard 0{as) contributions in PQCD. 

(b) 0{q:s) factorization into a wave function and a hard scattering amphtude. 

(c) Separation of hard and harder 0{as) contributions in an effective field 
theory, (d) 0{as) factorization into a "harder" function and a hard scattering 
amphtude. 

Fig. 2. (a) External W emission, (b) Internal W emission, (c) and (d) 
Nonfactorizable internal W emission. 
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(d) 
Fig. 1 




(d) 

Fig. 2 



